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In protocols of distributed quantum information processing, a network of bilateral entanglement 
is a key resource for efficient communication and computation. We propose a model, efficient both 
in finite and infinite Hilbert spaces, that performs entanglement distribution among the elements 
of a network without local control. In the establishment of entangled channels, our setup requires 
only the proper preparation of a single elected element. We suggest a setup of electromechanical 
systems to implement our proposal. 
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The role of entanglement in delocalized architectures 
of a device for quantum information processing (QIP) 
has been investigated under many aspects jj. Entangle- 
ment between distant sites of a distributed register is a 
fundamental requisite to optimize communication proto- 
cols and perform efficient quantum computation 2j. In 
this context, an entanglement distributor creates an en- 
tangled network of the elements of a register that, other- 
wise, have no direct reciprocal interaction. The efficiency 
of the distributor can be quantified by the number of el- 
ements which are entangled per single use of the distrib- 
utor or by the amount of entanglement shared by any 
two of them. Thus, the choice of the most appropriate 
design of the distributor is a problem-dependent issue 
with no general recipe. An interesting configuration for 
this problem is a star-shaped system in which an elected 
element interacts simultaneously with many other inde- 
pendent subsystems 53]- 

In this paper, we propose a model that acts as an ef- 
ficient entanglement distributor. An important feature 
of our proposal is that no local control on the dynamic 
of the participating systems is required once the interac- 
tions are set. We only need the pre-engineering of the net- 
work and a proper control of the interaction time. This 
is an advantage exploitable in those situations (frequent 
in solid-state physics) where single-element addressing is 
hard or impossible. The interaction we suggest acts on 
a multipartite bosonic network whose evolution can be 
tracked analytically both in the discrete and the contin- 
uous variable (CV) case. 

Despite our proposal being naturally described using 
the quantum optics language, we show that our model is 
general enough to find interesting applications in solid- 
state physics. We sketch a system of coupled electrome- 
chanical oscillators to embody our model. Similar setups 
have recently found applications in the entanglement- 
transmission problem [j| . 

The model - We consider N bosons (or modes) bj 
{j = 1, .., N) described by the annihilation (creation) op- 
erators bj (Sj) and an additional mode, labeled a, which 
we call the root. The interaction configuration is sketched 
in Fig. n (a) and consists of the resonant couplings of the 
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FIG. I: (a): The interaction configuration in Eq. Each 
edge represents an interaction, (b): Complete entanglement 
graph generated by Eq. 10. Solid and dashed edges represent 
entanglement. 



root to each bj. The satellite elements bj do not mutu- 
ally interact. In the interaction picture, we consider the 
Hamiltonian 
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+ h.c. (h^i) 
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with Qj real and time- independent couplings. For = 1, 
U{t) = e~^^''^ is similar to a beam-splitter (BS) super- 
imposing mode a to 61 . We analyze the characteristics of 
the many-body dynamics corresponding to A^ > 1. After 
a lengthy calculation based on Lie algebra we find that 
U{t) can be decomposed as 
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1 ■Rbjv-j^ftN-j + i6N-j(£Ar-J>0) 
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where = ELiS^ Rb, = e 



is a TT-phase shifter 



for mode bj, Sj = cos ^{Qj+i/'dj+i) and Bab{v,(p) — 
g[«(atbe'^-a6te-»'^)] ^^^^^^^ ^ Operator. This decom- 
position is extremely useful as it shows that the dynamic 
can be interpreted as the action of a setup of optical ele- 
ments on N + \ bosons. Eq. (jSJ describes how the root 
gains information from 6j's via the interaction Bb^a as 
well as the distribution of any information initially in a 
to the satellites. The form of Eq. reveals that, if bj's 
are all prepared in rotationally-invariant states (such as 
|0) or thermal states), the transformations prior to Bb^,a 
do not contribute to the entanglement dynamics j^. By 
properly setting the Sj , the evolution of the network can 
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be made equivalent to an array of BS's which sequentially 
superimpose a to bj modes. If a is in a superposition of 
|0) and a coherent state, 6j's being in the vacuum state, 
we generate an (iV+ l)-mode GHZ state useful for secret 
sharing 0. The entire Eq. ^ must be considered if we 
initially prepare one or more satellite modes in a coherent 
or a non-classical state. 

The model described by Eq. realizes various in- 
terference patterns in the equivalent all-optical setting 
allowing for different tasks. For instance, if Qj = Q 
(Vj), Hi describes an effective XF-coupling suitable for 
\ N phase-covariant cloning 0. As another ex- 
ample, let us take iV = 2 so that Eq. ^ reduces to 

ilir) = Rb,BbM£i,0)Bb,aiV2gT,-Tr/2)Rb,Bb,bA£i,0) 
with El = n/4. We assume that mode 6i is initially pre- 
pared in the single-excitation state |l)f,^, &2 and a being 
in the vacuum (the investigation can be generalized to 
the case of 6i being prepared in a coherent state). It is 
easily seen from our decomposition that at \/2Qt — n, 
the initial state is transferred to mode 62 with unit prob- 
ability (while the maximum probability of finding the 
initial state in a is only 1/2). This analysis shows that 
perfect quantum state transfer from bi to 62 can be per- 
formed through mode a. In fact, when ^/2Qt ~ mr 
(n = 0, 1, • • • ), it is interesting to note that our model 
is equivalent to a Mach-Zehnder interferometer with a 
mr phase-shift in the path of mode 62, which is obvious 
from our decomposition in Eq. 1^)1. 

Single excitation case - Consider a initially pre- 
pared in |1)^, bj's being in |0)j- The dynamics 
is captured by considering a fictitious collective mode 
of its annihilation operator c = ^j^j I'^n- Thus, 
C>M|10-0)„fe,..,„ - cos(z?Arr)|10),,-»sin(^^T)|01),, 
with 11)^ = ct \0..0)b,..b^ ~Y.j{Q3MN) |o..i..o),^..,^...,„. 
This state can be pictorially described by complete en- 
tanglement graphs as those shown in Fig.^(b). There, 
solid or dashed edges represent entanglement. 

In the basis {|00) , |01) , |10) , |ll)}f,,b^ , the reduced 
density matrix of the generic pair bi, bj i,j) reads 

/l-(Gf^+G2^) 0\ 

p . ^jN GiPfGjN 

GiNGjN GjjY 

\ 0/ 

where i < j and GjN = Qj shi{'dNT)/'d^. The entangle- 
ment of this mixed state can be quantified bv the concur- 
rence Cn — max{0, ai — 02 — 0:3 — 0:4} Here, a^'s 
are the square roots of the eigenvalues (in non-increasing 
order) of pij{(Ty (g) (Ty)p*j{cry (g) a-y) with p*^ the com- 
plex conjugate of pij and cry the y-Pauli matrix. We 
get Cn = max{0, 2GiArGjAr}. For later purposes, it is 
also useful to consider the entanglement measure based 
on negativity of partial transposition (NPT) 0. NPT is 
a necessary and sufficient condition for entanglement of 
any bipartite qubit state The corresponding entan- 



glement measure is defined as NPTm = max{0, — 2A } 
with A~ the negative eigenvalue of p^^j which is the par- 
tial transposition of pij with respect to bj. We find 
NPTn = max{0, [(1 - G^^ - G^^)^ + 4GfjvG|^]i/2 _ 
(1 — GjjY ~ ^^n)}- and NPTn arc optimized when 
i^nt = (2A:-|-l)7r/2 {k G Z). Using this condition as a con- 
straint in the Lagrange's method of indeterminate multi- 
pliers, we find that Cn and NPTn are maximized for the 
uniform set of couplings Gj — Q (Vj). In this case we get 

CN.v^a■. = 2/N and NPTN.macc = {[4+ (A^ - 2)2]i/2 _ 

(N - 2)}/N. 2/N is the upper bound for the bipartite 
entanglement in an iV-party system 0. Thus, Eq. ((U 
is optimal under the point of view of pairwise entangle- 
ment distribution. For equal Gj, the pij are all equal 
and we have |10..0)^(,^ ^ cos(??jvr) |1, 0..0)^j^ j,^ - 
ism{-dNT)\0,WN)ab b ■ have introduced the N- 

particle I^-state 11^^)," = N-'/^ |0..1..0),^..,^...,„ 
which is the state achieving GAr_maa;0|- Thus, the maxi- 
mum concurrence between any pair of bj's is found when 
the root is separable from the rest of the network. The 
corresponding graph is obtained by deleting the dashed 
edges in Fig. (b) , the satellite elements forming com- 
plete and permutation-invariant entanglement graphs. 
The system periodically evolves from a separable state 
to a configuration where the root is factorized from the 
rest of the network (which is in iWAr)^^ ^jv)- between, 
an {N -\- l)-partite entangled state is obtained. 

Recently, a configuration of many spin- 1/2 systems 
analogous to Eq. has been proposed Q. The one- 
excitation case we have considered allows for a compari- 
son between the two situations, both achieving CN,max — 
2/N . In our model the bosonic nature of the register al- 
lows for this result without local control on the satellite 
elements or the root. In ref. 0], on the other hand, this 
is obtained by using an additional magnetic interaction 
and through the measurement of the state of the root. 

In order to further characterize our entanglement dis- 
tributor, we compare IWOv)^^ bN ^'^ class of cluster 
states. These are known to be useful and genuine mul- 
tipartite entangled states |0, inequivalent to \Wn) for 
any N . While there are always proper local measure- 
ments on a subset of a cluster that allow for the deter- 
ministic extraction of a pure Bell state, this is not the 
case for a M^-state. However, the quantum correlations 
in a cluster are encoded in the system as a whole and any 
pairwise entanglement (obtained by tracing out the rest 
of the cluster) is zero. This is a drawback in those situ- 
ations where bipartite entanglement is required but the 
physical system is such that the realization of a measure- 
ment pattern is made difficult by the problems related 
to single-element addressing. Finally, the entanglement 
of IW^at) is persistent as N — \ projective measurements 
are required in order to disentangle the elements of the 
register. For the problem we address here, our analysis 
shows that Eq. J^l is a suitable and exploitable model. 
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CV case - Considering only the case of a single excita- 
tion in the root restricts the possibilities offered by the 
bosonic nature of our register. In ref. it is shown that 
a non-classical input is a fundamental pre-requisite for 
the entanglement of the outputs of a beam-splitter. The 
same is true in our case because of the analogy between a 
BS and Eq. On the other hand, necessary and suffi- 
cient conditions for the entanglement are known and en- 
tanglement can be quantitatively determined only for the 
class of two- mode CV Gaussian states 0| . In virtue 
of these considerations and because of the Gaussian- 
preserving nature of the linear operations in (0), only 
Gaussian states will be considered here. 

A powerful tool in the analysis of A^-mode CV systems 
is given by the variance matrix V, defined (after uni- 
tary displacements) as Vafs = {{xajXis}) {a, (3 — I, .., A). 
Here, x = {qi,pi, .., QNiPn)'^ is the vector of the quadra- 
tures. A Gaussian state is fully characterized by the 
knowledge of just the first and second moments of x and, 
in order to characterize the state of our A^ -|- 1 modes, we 
need to find the variance matrix of their joint state after 
U{t). In phase-space, the action of U{t) is such that 
V = T^~VT becomes the new variance matrix. Here, T 
is the 2{N + 1) x 2{N + 1) unitary matrix (found using 
Eq. ©) 



T = 



/cOs(t?ArT)l2 Ai<Ty A2<Ty 

Ai(Ty Diih Di2h 



AjMCTy 

DiNh 



(4) 

where I2 is the 2x2 identity matrix. An = —iGnN and 

Dnm = Snm + [^n^/m (cOs('!?Arr) - !)/'&%]■ Snrn deuOtCS 

the Kronecker symbol. 

For simplicity, we take Gj — Q (Vj), 6j's being in the 
vacuum state (variance matrix V^j^. b^ —®^=i^2bj)- The 
root is prepared in a squeezed state (squeezing parame- 
ter r) which is the most natural non-classical Gaussian 
state 0. The initial variance matrix of the system is 



6i..biv 



with V„ 



the variance matrix of a 



and (Tz is the z— Pauli matrix. By tracing all the modes 
but hi and bj we get 



-'N 



(5) 



Here, L^v = diag (n^r, mAr), Cat = diag(cAr, dA?) with 
UN = 1-1- Cat, niN = 1-1- (In, cn = —e^diq = (e*" — 
1) siT?{dNT)IN . No dependence on the indices i, j exists 
so that Eq. Q is the same for any pair. V'b^f,^. has a form 
which allows us to quantify the bipartite entanglement. 
Indeed, for a variance matrix as Eq. Q, the NPT entan- 
glement measure is given hy £n — maxjO , ((5i i52)~^ — 1} 
with 5i = TiAf — |cjv| and 82 = niN — |c?Af| |l3|. We have 



£n — max i 



2(1 - e-'')sin^(^?AfT) 
A - 2(1 - e-^)sm^ (i^NT) 



(6) 



which is plotted in Fig. 12 (a) against the effective cou- 
pling g = Qt. £n diminishes as A^ increases and, for 
fixed values of r, is maximized at i^nt = (2A: + l)7r/2 
{k E Z). In Fig. [21 (a) only £n>3 is shown as A = 2 
requires some comments. For this particular case, by 
generalizing the results of the analysis in refs. 0, Q|, 
we expect the evolved state of modes 61, 62 to be lo- 
cally equivalent to a two-mode squeezed vacuum. This 
result is crucially dependent on the fact that, from Eq. (O 
for A^ = 2 and i92T = 7r/2, the interaction between the 
satellite modes is an effective 50 : 50 BS. This allows 
us to decompose the variance matrix of the resulting 
two-mode state as 'V'bib2 = 0(^)S;,jh2(|)0(^). Here 
O(^) = (^jLi^bji^) with Sbj the single-mode squeez- 
ing transformation (which does not modify the entan- 
glement structure) and Sb^tj (|)_the variance matrix of 
a two- mode squeezed vacuum [l3|. The state is pure 
which implies the separability of a from bi -f 62. By 
studying the purity 'Pbib2 = [det V'{,^b2]~ = ^sji we find 
that its period is one- half the period of £2- That is, the 
bi + 62 state is pure not only when £2 is maximum (at 
Todd = (2fc-|-l)7^/(2^/25)) but also at Teven = kir/iV^G), 
which corresponds to £2 = 0. By using the biseparability 
condition of a boson from a group of A^ others |l(il |. we 
have also checked that at Teven no entanglement is found 
between a and bi + b2- The state is fully separable. 

By enlarging the network to A^ > 3, we notice that the 
first interaction in Eq. ^ between two satellite modes is 
a BS of its Sn-1 — 1 /VN, which is no longer a 50 : 50 BS. 
This stops the possibility of getting a state locally equiva- 
lent to a two-mode squeezed state ^3 ^^'^ structure 
of the multi-mode entangled state becomes much more 
complicated than the simple case of A^ = 2. In particu- 
lar, the state of any pair (bi, bj) is pure just at geven but 
no longer when £jq is maximum. Thus, quantum correla- 
tions are shared between bj's at Todd = {2k+l)7r / {2^/NG) 
but not between the root and them. The entanglement 
configuration alternates between a fully separable state 
and a many-body entangled state of just 6j's, passing by 
a configuration in which entanglement is shared with a. 
The picture given by the graphs in Fig.n(b) is still valid. 
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FIG. 2: (a): £m against the dimensionless time g = Qt for 
iV = 3 (solid line), iV = 4 (dashed line) and = 5 (dot- 
dashed). The squeezing of the initial root state is r = 0.8. 
(b); Relative entanglement differences Ai (if) and Acv (♦) 
against A''. 
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We now look at the effect of increasing N on the prop- 
erties on the entanglement distribution. We consider 
the quantities Ai = 1 - {NPTpf+i,jnax/NPTN,jnax) and 
^cv = 1 - {£N+i,max/£N,max) which measure the rel- 
ative loss in pairwise entanglement if the network is en- 
larged by one element. Fig. 12 (b) shows that at a fixed 
r, Ai and Acv decrease with iV (3 < < 20). The dis- 
tribution process is only weakly affected and the entan- 
glement is still spread through the network. In passing, 
it is interesting to stress the qualitative robustness of the 
distributed entanglement in the CV case as compared to 
the discrete one, an issue which, in a different context, 
has also been noticed in [l8j |. 

Possible setups - We briefly mention that, to embody 
Eq. 1^, we can use the interaction of a linearly polar- 
ized optical bus with N ensembles of cold atoms (con- 
fined in vapor cells), providing the Hamiltonian Hie — 

i^PphJ^i^iPei (^^ is a couphng rate). Here, Pph (pei) is 
the momentum operator of the bus {i*^ atomic ensemble) 
whose wavelength is assumed to be much larger than the 
dimensions of the ensembles and their separations [l5| . 
Hie holds within the Stokes-vector formalism for the bus 
and the Holstein-Primakoff transformation mapping col- 
lective states of an ensemble into a fictitious boson. By 
discarding rapidly-oscillating terms. Hie Hj. 

Stimulating opportunities come from micro and nano- 
electromechanical systems (MEMS and NEMS), i.e. elec- 
trically controlled mechanical oscillators (or cantilevers) 
whose dimensions are in the range from 10~^ to 10^^ m. 
Doubly clamped cantilevers with fundamental mode fre- 
quency in the range of [10^ — 10^] Hz have been fab- 
ricated and mutually coupled 0. They are useful to 
study Heisenberg-limited measurements and entan- 
glement 0, 01 . There are theoretical proposals for 
ground-cooling and squeezing of NEMS mode . The 
preparation of phonon-number states and the tomogra- 
phy of a vibrational mode have also been addressed |22l | . 

We consider N classical oscillators coupled via spring- 
forces to a central one, the analogue of our root. Within 
Hooke's law, the energy of the system is Ti = {uj/2)[q^ + 
pI) + (w/2) I?' +Pj + ^0 ilj - 9a)'] , where the JC^ 's are 
the coupling factors, {qj, pj) are proper canonical vari- 
ables and u! is the frequency of the oscillators (equal for 
all). Each ICj is controlled via voltage biases between the 
cantilevers. Each bias creates a potential that changes 
with the capacitance between two oscillators. Eq. is 
then found in a second-quantization picture and within 
the rotating wave approximation (used for ICi ~ O.lo;). 
The oscillators can be built via photolitography of gold 
on silicon substrates [13|. In our case, planar grids of a 
few cantilevers face each other in pairs, surrounding the 
root. The coupling of the cantilevers to the phononic 
modes of the substrate is the main source of decoher- 
ence. However, oscillators with quality factors Q ~ 10* 
and Lu ~ 10 MHz (coherence times ~ 1 msec) allow 



now for a good number of coherent operations. The re- 
construction of VJj.j. is challenging here. However, a 
single-electron-transistor (SET) capacitively coupled to 
the cantilevers can be used |21|. Exploiting the changes 
of the coupling capacitances (which depend on the in- 
stantaneous position of the oscillators), a SET acts as 
a displacement-to-current transducer with displacement 
sensitivity ~ 10~^^ m/yj Hz. Stroboscopic techniques to 
infer VJ,.^^ could then be used 0. 

Remarks - We have characterized a many-body interac- 
tion that, through just global interaction with a seeding 
system, distributes entanglement in a network of local 
processors. The dynamics is described by linear opera- 
tions and the model is flexible enough to allow for differ- 
ent interference patterns by pre-engineering the couplings 
and the initial state. We have shown how symmetric bi- 
partite entangled states are generated both in the discrete 
and CV case. To embody our model, we have described a 
setup of coupled cantilevers that offers nice perspectives 
in the study of entanglement distributors for QIP. 
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